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The aim of this note is to show that the results of the recent experiment 
measuring EPR-type flavour entanglement in Y (4S 1 ) — > B°B° decays^ are 
compatible with local realism. It is known that a Bell inequality test cannot 
be performecP but this does not prove that the experiment is compatible 
with local realism. I shall prove the compatibility by exhibiting a local hidden 
variables (LHV) model which reproduces the quantum prediction (and agrees 
with the obtained results within experimental errors). The time-dependent 
rate for decay into two flavour-specific states are^ 

Ri(At) = exp(-At/r) (l + (-1)* cos (AmA*)) , At=\t 2 -h\, (1) 

where Am is the mass difference between the two B° — B° mass eigenstates, 
i — 1 corresponds to the decays B°B° or B°B° and % = 2 to the decays B°B° 
or B°B°. Actually Ri are the probability densities of decay at time t 2 of 
the second particle (say the one going to the right) conditional to the decay 
of the first (say going to the left) at time t\, both t\ and t 2 being proper 
times of the corresponding particles. For our purposes it is more convenient 
to consider the joint probability densities, rki{t\,t 2 ) for decay of the first 
particle at time t\ and second at time t 2 , where k — 1 (I — 1) means that the 
first (second) particle decays as B° and k = 2 (Z = 2) means that the first 
(second) particle decays as B°. 

According to Bell's definition of LHV modelj^ appropriate for our case, 
we should attach hidden variables Ai and A2 to the first and second particles, 
respectively, and define probability densities p, Pk, Qi such that 

r u (ti,t 2 ) = / p (Ai, A 2 ) P k (A 1? h) Q l (A 2 , t 2 ) d\ x d\ 2 . (2) 



The function p, giving the initial distribution of the hidden variables in an 
ensemble of Y (45) decays, should be positive and normalized, that is 

P (Ai, A 2 ) > 0, f p (A 1; A 2 ) dX l( iX 2 = 1. (3) 



The functions Pk (\i,tx) represent the probability density that a particle with 
label Ai decays at time t\ as a B° (^B ^ if k = 1(2) and similar for Qi. Thus 

these functions should be positive and, as all B° or B° particles decay sooner 
or later, they should be normalized for any {Ai, A 2 }, that is 

/»oo 2 /»oo 2 

Pk,Qi>0, dt 1 J2P k (\ 1 ,t 1 ) = l, dt 2 J2Qi(^2,t 2 ) = l. (4) 
Jo k=i Jo 1=1 
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Any choice of functions {p,Pk,Qi} fulfilling eqs.(TS]) to (jlj) provides a LHV 
model predicting the joint probabiliy densities of decay r ki (ti,t 2 ) ■ 

I propose the following. For the initial distribution of hidden variables 

P (Ai, A 2 ) = At ^ X2) S (Ai - A 2 ) , A 1; A 2 G [0, 2tt], (5) 

where 5 () is Dirac's delta, and the functions A" (A 2 ) will be defined below, 
after eqs.(JHJ) ? where the normalization of p(X\, A 2 ) will be proved. For the 
probabilities of decay 

1 / t \ °° 7T 

PkiKh) = - exp (--M 0(ti)^e (~- |Ai + (2n-fc)7r- Amti|) , 

(6) 

<5/ (A 2 , t 2 ) = iV (A 2 ) exp f-^j (t 2 ) [cos (A 2 - (/ + 1)tt - Amt 2 )] + , (7) 

where (t) = 1 (0) if i > (t < 0) and means putting if x < 0. Thus 
all four functions are decaying exponentials modulated by periodic funcions 
which oscillate with period 27r/Am. Physically this means that each particle 
"lives" as a B° during a time interval of duration n/Am, then becomes a 
B° during another time interval n/Am, and so on, until it decays. The 
particles are always anticorrelated in the sense that, at equal proper times, 
one of them is B° and the other one is B°. 

From eqs.(|Sl) and © it is easy to see that the total probability densi- 
ties (i. e. independently of flavour) for the decay of particles 1 and 2, are 
respectively 

I>(Ai,ti) = 7 ex p(-|)' 
fc=i ' ^ ' 

^Qi (A 2 ,t 2 ) = N(\ 2 ) exp |cos(A 2 - Amt 2 )|, 

i=i V r / 



(8) 



where ti,t 2 > 0. We see that the decay of the first particle is given by 
a standard exponential, but the decay law of the second particle is more 
involved. The functions A^ (A 2 ) are chosen so that the normalization eq.fll]) 
holds true. It is not necessary to calculate explicitly the functions N(X 2 ), 
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which are rather involved, but I derive an important property, namely 

/ — ——d\ = dX exp I — I |cos(A — Amt) \ dt 

Jo N{\) J J V t) 

= jf exp y--J dt \cos(\- Amt)\d\ = 4t. (9) 

This relation proves that the distribution p(Ai, A 2 ) , eq. (JSJ) , is indeed nor- 
malized. 

In order to get Tki (ti,h) we should insert eqs.Q and (GO) in eq. (J2|) and 
perform integrals which are straightforward. Introducing the new variable 

x = Ai + 2mr — kn — Amti, (10) 

and performing the integral in A2, using Dirac's delta, we get 



( . . s 1 ( h + 1 2 

rki{ti,t 2 ) = ^2" ex P 



Iki = dx [cos (x + (A; — / — 1) tt + s)] + , s = Am(ti — 

J-tt/2 

where we have taken into account that only one term of the sum in n may 
contribute, depending on the values of t\ and t 2 , and we have removed the 
irrelevant term 2nir in the argument of the cosinus function. It is easy to see 
that the functions Iki are periodic in the variable s with period 2tt. Thus it 
is enough to consider the interval s G [0, 2ti] . Thus in the particular cases 
k = l = lork = l = 2 the integral (TTTT) becomes, for s G [0, tt] 



hi (ti, t 2 ) = I 22 h) = / [cosx] + <ix = / cosxda; = 1 — cos s, 

Jw/2-s J-k/2-s 

(12) 



and for s G [7r, 27r] 

/•3tt/2—s p3tt/2—s 

hi (h, ^2) = I22 (ti, t2j = \ [cosx] + dx= / cosxdx = 1 — coss. 

Jtt/2-s J-tt/2 

(13) 

Similarly we get, for any s = Am (ti —t 2 ), 

h2{h,t 2 ) = hi{h,t2) = 1 + coss, (14) 
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Finally we obtain 



r M (h, t 2 ) = -L exp \Jl±±*)j [1 - (-l) l - k cos (Am At)], At = \t t - t 2 \ . 

(15) 

Hence we may get eq.([T]) via the equality which defines the conditional prob- 
ability reported in the commented paper^ in terms of the joint probability, 
namely 

Ri = rexp(2tj/r)r fe/ , j = \k - l\ + 1, 

where i = 1(2) if particle 1(2) is the one decaying first. This proves that our 
LHV model's prediction agrees with the quantum one for the said experi- 
ment. 

The model may be interpreted physically saying the either particle pro- 
duced in the decay of the Y(4S) oscillates between the two flavour states in 
such a way that the flavours or the two particles in a pair are opposite at 
equal proper times. The model looks somewhat contrived due to the lack 
of symmetry, in the sense that the functions P^ are quite different from the 
functions Qi. A more symmetrical model may be obtained assuming that 
the assignement of the functions Pk and Qi to the particles in a pair is at 
random. In any case our purpose was only to show the compatibility of the 
experiment with local realism, and not to make a physically plausible model. 

I acknowledge useful comments by Albert Bramon and Alberto Ruiz. 
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